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Abstract 

We report direct observation of the evolution of local temperature inhomogeneity and the result- 
ing atomic migration in a metal film (Ag on Si) stressed by a current by using a Scanning Thermal 
Microscope that allows simultaneous temperature mapping and topography imaging. The experi- 
mental observation is analyzed using a model based simulation. The experimental observation and 
the simulation show that due to current stressing the temperature of the film becomes significantly 
inhomogeneous over time (with local temperature deviating strongly from the mean). This creates 
local stress as well as local temperature gradient that lead to mass migration in addition to the 
electromigration. We show that the local temperature inhomogeneity serves as one of the main 
agents for local atomic migration which leads to change in film microstructure. The migration leads 
to damage and eventual failure as simultaneously monitored by in-situ resistance measurement. 
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I. INTRODUCTION 



In presence of large electric current through a metal film, ions in the film are acted upon 
by a number of forces (in addition to the direct electrostatic force) that lead to migration 
of ions . These forces can be due to momentum transfer from electrons to the ions leading 
to electromigration (EM) or due to stress (stress migration) or stress resulting from thermal 
gradient , |6j|. Whatsoever may be the cause of this mass migration, it leads to 

change in film microstructure and formation of voids or hillocks leading to complete damage 
when the current density is very high. In this paper we report an important experimental 
result which clearly shows that the initiation of mass migration in a current carrying metal- 
lic film is associated with development of local temperature inhomogeneity (spatial) which 
acts as the sites for nucleation of atomic migration and further microstructural evolution 
depends substantially on the evolution of spatial temperature inhomogeneity during the en- 
tire damage process. This shows that in addition of the average temperature of the film, the 
local temperature (which is different from the average) plays an important role. We did this 
by imaging the temperature variation of the surface of a current carrying metal film using 
a scanning thermal microscope (SThM) which allows us also to generate a simultaneous 
topography map in combination of in-situ measurement of the resistance of the film. It has 
already been demonstrated that SThM can be used as a powerful thermal imaging tool to 
map the temperature distribution due to joule heating effect in metal lines carrying high 
current [zj]. The dependence of atomic migration through metal films on the equilibrium 
temperature of the film and the joule heating effect causing change in average film tem- 
perature have been investigated before j^l, 9]. However, to understand the dependence of 
the damage process on thermal distribution, it is very important to study the correlation 
between evolution of microstructure and the spatial inhomogeneity of temperature as a func- 
tion of time. In particular, development of local temperature gradient (VT) and its impact 
on the mass migration are important factors and it will be important to follow the evolution 
of these factors as the damage process evolves with time. To our knowledge, this impor- 
tant issue has not been addressed to in any earlier study. Our experiment, being based on 
simultaneous topographic imaging and spatially resolved temperature mapping along with 
in-situ resistance measurement with the progress of the da mag e process, directly addresses 
these issues. Imaging techniques like Infrared Microscopy lOj, Laser- reflection thermome- 
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111 ]. SThM [7] and Scanning Joule expansion microscopy [12] have been used in the past 



to image the thermal inhomogeneity of metal interconnects with different degrees of spatial 
resolution. However, the evolution of the inhomogeneity, local temperature gradient and 
its direct correlation to the mass migration as the damage process evolves (established by 
in-situ monitoring of the resistance) have not been done before. In addition to the direct 
experimental observation we also carry out a simulation based on a model to analyze the 
effect of the local temperature inhomogeneity on the mass migration process and thus gain 
understanding of the observed results. 



II. EXPERIMENTAL TECHNIQUES 

Experiments were carried out in Ag films (thickness ~ 0.160/im) as test samples (dimen- 
sion pa 20 fim square). The films were prepared by thermal evaporation on Si (held at 175°C) 
pressure of 10 8 mbar and post-growth annealing at 250°C for 6 hours in the same 
vacuum. The XRD of the films show that the films were mostly oriented in (111) direction 
and well textured. The room temperature resistivity (p3oox) and residual resistivity ratio 
(P300K / P4.2k) for the films were typically around 1.65/xf2cm and ~ 8 respectively. 

FIG [1] shows the schematic for in-situ resistance and SThM as well as topography 
measurements. The SThM and topography measurement was based on a contact mode 
Atomic Force Microscopy (AFM) technique with an additional temperature measurement 
unit (SThM module) and using the cantilever of AFM itself as a thermometer. We used a 
commercial AFM from Veeco Metrology Group ( Model: Autoprobe CP Research) and a 
cantilever with a built-in resistive thermal element made of Pt/10% Rh alloy which acts as 
a probe for scanning thermal microscope |l3||. A mirror, cemented on the cantilever, acted 
as the reflector for the laser light to measure the deflection of the cantilever. To measure 
the temperature of the sample at the point where the tip of the probe touches it, the probe 
was connected as an arm of a balanced Wheatstone bridge powered by a constant current 
from SThM module. The error voltage, V out (see FIG [1]), from the bridge varies linearly 
with local temperature (T) of the point where the tip touches the film and is given as 

V out = Vo + taR C(T-T ) (1) 

where, Vq = bridge voltage at the reference temperature, i = current applied to the probe, 
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a = temperature coefficient of resistance of the probe material, Ro = resistance of the 
probe at a reference temperature To and C = calibration gain of the SThM module. Before 
carrying out measurement on the actual sample, the SThM probe was calibrated against 
V out . This was done by using a PtlOO film thermometer as a sample which was placed on a 
heater and recording V out for different reference temperatures as measured by the PtlOO film 
thermometer which is serving as the sample. From a fit through these calibration points the 
absolute temperature could be calculated according to the equation, 

T = (T - -V Q ) + —Vout (2) 
m m 

for any value of V ou t- From the calibration curve we find the slope m = 8.3mV/K. The 
spatial resolution of the thermal (local T) images was limited by the finite size of the tip 
to ~ 200nm. The test sample size was kept in the range of ~ 20/xm due to this spatial 
resolution . We applied a direct stressing current (dc) of density Jd c = 2 x 10 7 A/ cm 2 across 
the current leads and the voltage drop was recorded to measure the resistance. The damage 
process was monitored along with imaging by directly recording the resistance (R) at regular 
time interval automatically until the film got damaged. With this experimental setup we 
could investigate systematically the evolution of the damage process with in-situ AFM and 
SThM. 



III. RESULT 



In the figure [2] we plot the evolution of the resistance of a film as a function of time with 
the progress of the damage process. Resistance is plotted in logarithm scale to enhance the 
initial lower resistance regions. It was observed that initially the resistance increased very 
slowly for first 105 hours after which (marked as ti in the plot) it showed rapid increase. 
Before reaching the point t\, incubation of resistance followed R = Roe Plt kind of trend, 
where, Rq is the initial resistance of the film and p\ ~ 2 x 10 _3 /i _1 is the slope of the 
logR — t plot. After this point also increase in resistance followed a similar trend but with a 
much larger slope P2 ~ 3.6 x 10 _2 /i _1 , showing a rapid increase. This trend was seen to be 
maintained till around 120 hours (marked as t 2 in the plot) after which also resistance kept 
increasing at a similar average rate. In the late stage of damage process, as shown in the 
FIG [21 there are steps in resistance changes. We generally take the SThM and topography 
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images in these regions because the resistance stays constant over a time period. The steps 
are marked as to, t\ ti in the figure. 

In the FIG [31 three representative topography, SThM images and temperature profiles 
are shown. The data obtained at different stages of current stressing process, are shown in 
three different rows. The progress of the damage process can also be seen from the change in 
R of the film with the time of current stressing (FIG [2]). Resistance change in the FIG [2] and 
the AFM based imaging were taken simultaneously. The first figure in the row (a) shows 
the topography of the whole film including contact pad, while in all the other images only 
the part of film, marked by a square, is shown to elaborate the changes that occured as 
a result of current stressing. Row (a) corresponds to the film at unstressed condition (at 
to)- The line-scan is a typical example of temperature inhomogeneity in an unstressed film 
and shows the noise limited resolution. As the damage process progresses one sees there 
appearance of relatively hot areas surrounded by relatively cold contours as indicated by 
arrows in the SThM image in row (b), which was taken after 146 hours (at t$) of current 
stressing. A typical line-scan on SThM image shows that the spatial temperature fluctuation 
has gone up. Simultaneous contact mode AFM mapping shows that these local hot spots 
act as the regions where the Alms's microstructure has changed considerably. This becomes 
more visible (see row (c)) in later stage of the damage process where there were sufficient 
changes in microstructure and the film becomes rough (with formation of hillocks and voids) 
at these positions and a more non-uniform temperature profile develops. The row (c) images 
were taken after 163 hours (at £5) of current stressing. As the damage process progresses, 
the temperature over the film become more non-uniform and a larger part of the film gets 
damaged. [Note: The line scans, shown in the FIG\^at different stages of current stressing, 
were taken at the same sites to show gradual evolution of temperature inhomogeneity. 

The result of in-situ SThM measurement qualitatively shows that the progress of damage 
process is accompanied by the increase in surface roughness, enhanced local temperature 
variation and also an enhanced average temperature. A quantitative evaluation of these 
observations can be obtained from the development of these parameters with the progress of 
the damage process. Correlations among these parameters will lead to a more clear picture. 
To achieve this, we calculate rms roughness of the film (o~(h)), the average temperature 
(< T >) and the relative variances in temperature (< (AT) 2 > / < T 2 >) for an arbitrarily 
chosen part of the film at different stage of the damage process. The averages are taken over 
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the whole surface of the film. We have marked such an area with a square in the topography 
image in the row (a)of the FIG [31 The dimension of this area is 12/xm x 12/im. All the 
further calculations presented here were done for this area. 

To calculate the rms roughness of the film we find out the hight profiles along some lines 
over the selected part of the film. The standard deviation of the height profile gives the 
rms roughness along a line. We draw 60 numbers of such line in horizontal direction and 
60 lines in the vertical direction. Thus each line is separated by w 200nm, a value close to 
the resolution of the SThM tip. For each stage of damage we calculated the rms roughness 
averaged for all the lines. This rms roughness (averaged for the area) for the stages 

marked as to, t\, t 2 ^7 in the FIG H] are plotted in the FIG0Ja) as a function of time along 

with the resistances of the film at the respective stages. From this plot increase in cr(h) 
with the time of current stressing can be clearly seen. Since the roughness of surface has 
resulted from the mass flux in the film due to the migration process, the quantities cr(h) at 
any stage is a measure of the mass migration of the film that occured at the corresponding 
stage due to the current stressing. Initially cr(h) increased very slowly, but became very 
rapid when the resistance change is also rapid. The simulation, described later on, shows 
that this happens when there is significant mass migration. (In the example of the film 
shown here this happens after 146 hours of stressing). Towards the end of the experiments 
the cr(h) increased by more than one order of magnitude. 

To study the evolution of average temperature with the progress of the damage process, 
we calculate mean temperature (< T >) for the same part of film in the same way as was 
done for calculation of a(h). Here we take the thermal map of this part of film and find 
out the temperature profile along 60 numbers of horizontal lines and 60 numbers of vertical 
lines. Some typical temperature line profile are already shown in the FIG [3j The mean 
temperatures (< T >) for various stages of damage, calculated from the mean temperature 
for each line and averaged over the all lines give average temperature of this film area at 
corresponding stages of damage. The evolution of < T > is plotted in the FIG BJb) as a 
function of stressing time. With the advances of the stressing time < T > increases, initially 
little slowly but rapidly after 120 hours of stressing. Relative variances of temperature (< 
(AT) 2 > / < T 2 >) for any stage was obtained from calculating the variance of temperature 
for each line and averaging over all the lines and dividing the resultant quantity by the 
average of the squares of temperature along all the lines. The quantity < (AT) 2 > / < 
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T 2 > represents spatial inhomogeneity of temperature over the selected area irrespective of 
the background temperature. Change in relative variances in temperature is shown in the 
FIG^a) as function of time. Corresponding resistance values are also shown there. Almost, 
throughout the stressing time, all these quantities were seen to increase with advances of 
the damage. 

IV. DISCUSSION AND SIMULATION 
A. General discussions 

The results of the experiment show that local temperature variation that is created on 
current stressing is an important physical parameter in the current stressed film and it plays 
an important role in th eventual failure of the film. The increase in average temperature 
of the film is due to the joule heating effect. Due to microstructural change in the film 
caused by atomic migration during the damage process, the local resistance changes and 
this leads to spatial variation of the current density. This will lead to local variation in heat 
dissipation which will give rise to the observed spatial variation of local temperature in the 
film. This particular issue that local resistance inhomogeneity leads to inhomogeneity in 
local heat dissipation and hence the local temperature variation is also borne out by the 
model simulation discussed below. With the advances of the stressing time, increased surface 
roughness of the film with introduction and development of voids or hillocks in the film, such 
temperature distribution is expected to show stronger spatial variation, irrespective of the 
background average temperature. This is what was observed in this experiment. Comparing 
evolution of < AT 2 > / < T 2 > with that of cr(h), it could be seen that for most of the 
stressing time, relative spatial variation of temperature increased along with the increased 
surface roughness which became very rapid towards the end of the damage process with 
increase in value by more than two orders of magnitude. After the damage commences the 
local mass migration that results from the local temperature variation leads to enhanced 
average temperature (because larger heat dissipation resulting from larger R) as well as 
enhanced local gradient leading to more mass migration. Eventually there is a "run away" 
situation in all the parameters and the film is completely damaged. 

In a current stressed film, temperature has a great influence on the mass migration and the 
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local atomic flux. Particularly, atomic migration resulting from local temperature gradient 
or thermomigration becomes one of main components of the net mass migration, while the 
other components are electromigration and stress- migration. The component of mass flow 

-Eg 

arising from thermal gradient scales with i VTe k P T where E n is the activation energy for 

n 

diffusion [14|. The temperature map obtained from the SThM images at different stages of 
the damage process allow us to calculate the local temperature gradient. We obtained the 
magnitude of the temperature gradient (| VT |) by averaging over the same regions that 
we used to calculate rms roughness. To calculate | VT | we used same temperature profiles 
for all the lines used as described above. We obtained W- and |r from the horizontal lines 

ox ay 

and vertical lines respectively for each of their cross-point and then we calculated the local 
value of | VT | at that point by taking square root of (f^-) 2 + (f^) 2 - When this quantity 
was averaged over all cross-points, | VT | gave the magnitude of the temperature gradient 
averaged for the selected region of the film. In the FIG [5](b) we have plotted | VT | as 
function of time. To show development of surface non-uniformity with the evolution of 
temperature gradient, we have plotted cr(h) also in the same figure. It can be clearly seen 
from this plot that the local gradient, | VT | increases monotonically with time and that 
the rms roughness follows the same trend. Thus, the local temperature distribution as well 
as the local gradient are related and both have significant effect on the migration process. 
This particular observation that local temperature gradient controls in a significant way the 
direction of local direction of mass migration is an important out come of our experiment. 
The model described below provides support to this explanation. 

The effect of local temperature distribution on the mass flow can be mediated via some 
other related underlying physical processes. Particularly local variation in temperature 
produces a mechanical stress field in the film due to mismatch of thermal coefficients of 
expansion of the film material (Ag) and the substrate (Si). The stress field generated thus 
also acts as an additional strong driving force for atomic migration. Though it is not possible 
to single out the measure of damage resulting from the temperature gradient field alone, the 
results of in-situ SThM have shown clearly that local temperature variation has strong 
influence on the local atomic mass flux. 
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B. Model of the process and simulation 



The damage process of current stressed metal film is highly complicated process where 
various phenomena dependent of number of parameters like current density, operating tem- 
perature, thermo-mechanical properties of the film and other materials in physically contact 
with it. It is also important to note that the film cannot be considered as smooth. It has 
topological features and also local variation in resistance. It is thus not possible to obtain 
an accurate analytical result on the evolution of damage process . However, with a good 
model the damage process can be simulated. As a part of our present investigation, we 
used a realistic model and carried out a simulation in order to gain an understanding of 
the microstructural evolution as well as the local temperature evolution of the film as the 
resistance chang GS £lS cl result of the current stressing. We will see that the model, though 
simple captures most of the features observed in the experiment. Below we describe the 
model briefly and present the results. The detailed model and the computational methods 
are beyond the scope of this paper and have been presented separately lo) ]. 

In this simulation we used the initial topography of the film as revealed by the AFM 
image as the input. This was done by taking AFM topography image of pristine film. Then 
the film was modeled as a resistance network (see the FIG [6]) and the film/substrate system 
was modeled as finite element (FIG [7]). The value of the each local resistance element was 
obtained from the height profile of the AFM image. We made the assumption that on the 
average the resistivity is homogeneous across the film but the local resistance variation arises 
from the local height variation (that then controls the cross-section of the local resistance 
element) which is captured by the actual AFM image. The equivalent resistance of the film 
was obtained by solving the network. This equivalent resistance was then compared with 
the experimentally measured resistance of the film and a scaling factor was obtained. The 
scaling factor was then used to scale the values of all the resistance elements of the network. 
In this way we got a more realistic value of the local resistance elements that comprise the 
network and that is compatible with the resistance of the actual film. 

When the network is biased with a constant current source, the current through the 
each resistance element was calculated. Each resistance element also behaved as a source 
of local heat input due to joule heating. The finite element modeling of the film was done 
to calculate correctly the local temperature at any point of the film. The height of a finite 
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element for the film at any location was taken according to the film's height at that location 
whereas the height of all the elements for the substrate were kept fixed throughout the 
simulation. Using finite element analysis (FEA, described briefly in the appendix) the local 
temperatures for each node of all elements were calculated for the film and the junction 
between the film and the substrate. The corresponding stress due to mismatch of thermal 
expansion coefficients (TEC) of the film and the substrate were calculated. Using the just 
obtained local temperature and stress, net mass flux arising from combined contribution 
of 'electromigration', 'stress-migration' and 'thermo-migration' was calculated. For this we 
used the formula: 

Ja = ^f( Z *aepJ - - Jl-VT)exp[-E a /kpT] (3) 

where, J a is the atomic flux, D Q is the prefactor of diffusivity, C is atomic concentration, Z* 
is effective valence of the material, p is resistivity of the film, J is current density, Q is atomic 
volume, cr is the stress, Q* is heat of transport, Na is the Avogadro's number and rest of 
the symbles have usual meaning. For Ag film, we took values of D = 6.4 x 10~ 7 m 2 / s Ji| 
for self surface- diffusion, Z* = —21 p = 1.65 x 10~ 8 Q.m (experimentally obtained), 
Q = 1.7 x 10" 29 m 3 , Q* = 0.5 eV and N A = 6.022 x 10 23 . We used J = 2 x 10 7 A/cm 2 and 
kept E a as adjustable parameter. 

The resulting mass migration was then used to recalculate the new height profile of the 
film. The new local height profile becomes the input for the next simulation steps. At 
the end of each simulation step we can obtain the new resistance of the film, the average 
temperature, the local height variation (which we expressed as the rms roughness) and the 
local temperature variation that we expressed as the rms temperature fluctuation averaged 
over the film surface. These are also the experimentally determined observable. This process 
was followed in every simulation step. Thus, in our simulation we could see evolution of 
microstructure and local variation in temperature along with evolution of the film resistance 
as a function of time (which is the number of simulation steps) as the damage process progress 
result of the current stressing;. 

The results of the simulation (plotted as number of simulation steps) are shown in the 
FIG M and FIG |9j The results of simulation can be compared with the experimental results 
shown in FIG H] and FIG [51 In the FIG E^a) and E(b) we show the evolution of rms rough- 
ness (u(h)) and the average temperature (< T >) alongwith the evolution of resistance as 
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obtained from the simulation. In the FIG [J^a) we show the evolution of resistance of the 
film along with the change in relative variances of temperature. The evolution of surface 
roughness and that of the temperature gradient can be seen from the FIG[9](b). The simu- 
lation showed that the surface roughness, mean temperature, temperature gradient and the 
resistance increase slowly with time initially , but at very rapid rate after some instance 
when the resistance starts to change at a rapid rate. This result is in conformity with what 
has been observed experimentally.The simulation not only reproduces the experimental ob- 
servation qualitatively it also has quantitative validity. As an indicator we can compare 
the numbers when as result of the migration the resistance R of the film reaches 1Q. In 
experimental observation the parameter surface roughness <j(h) reaches ~ 0.017 fim. The 
simulation shows a value of ~ 0.023 /im. The average temperature < T > at that time is 
~ 339 K. The corresponding number from the simulation is ~ 321 K. The average local 
temperature gradient | VT | at that instance obtained from the experiment and from the 
simulation are 0.95 K{im~ x and 1.13 K^im~ l respectively. The Relative variances of tem- 
perature (< (AT) 2 > / < T 2 >) for that stage from the experiment is 2 x 10~ 6 and from the 
simulation it is 4 x 10~ 6 . The two local thermal parameters obtained from simulation though 
similar to the experimental result are however not the same. This reflects the complexity in 
modeling the underlying local heat transport process. Nevertheless we find that the simula- 
tion reflects in a significant way the process of current stressing. This is encouraging because 
there are no adjustable parameters in the simulation except the activation energy E a . 

The simulation establishes the basic understanding that due to increase in the local 
temperature gradient, enhanced mass migration was observed resulting in an increase of 
surface roughness also. The simulation shows that at this instance there are substantial 
mass migration in the film. The simulation shows that this increases the resistance values 
of the local elements leading to an enhancement of the average temperature are controlled 
largely by the local temperature gradient as well as the local temperature inhomogeneity as 
expressed in the local rms temperature fluctuations (in addition to the average temperature 
T). 

The simulation shows that the temperature enters in two ways. First, since the migration 
is an activated process, enhanced average temperature accelerates the process. Second, an 
enhanced local gradient of T also enhances the atomic flux. 
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V. CONCLUSION 



To summerize, in this investigation, we have carried out a systematic study on time 
evolution of local temperature variation as well as the local topology variation of a thin 
metal film subjected to high current stressing using AFM and SThM. We also measured 
in-situ the resistance change to monitor the progress of damage of the film under current 
stressing. We found that although the film has uniform temperature initially at some point 
of time there appear regions with relatively large temperature inhomogeneity with hotter 
zones surrounded by colder contours. The migration process nucleates in these regions and 
significant atomic migration results here as the time progresses. As the time progresses, new 
hotter and colder zones with higher local temperature difference appear. Our experiment 
demonstrates a close correlation among the atomic migration, the variation of local temper- 
ature and microstructural variation like surface roughness. The experimental observation 
was analyzed using a model based simulation . The results of the simulation show that the 
local heating controls the stress as well as the thermal gradient (| VT |). The migration 
due to local stress and the thermal migration add upto the electromigration induced mass 
migration and create the total the mass migration. This leads to damage of the film. Our 
investigation shows the predominant role of not only the average temperature growth but 
also of the local spatial temperature variations that determines the nature of net mass flow 
into or out of a region of the film under current stressing leading to its damage and eventual 
failure. 
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APPENDIX: A 



Finite element analysis was aimed to calculate temperature distribution throughout the 
film and even in the substrate at every simulation step by solving the heatflow equation 

dT 

p m s— = V.(kVT) + Q (A.l) 

where, p m and s are the mass density and the specific heat of the material respectively, 
T is the temperature, k is the thermal conductivity of the material and Q consists of the 
rate of heat generation per unit volume (as a result of Joule heating) and radiative heat 
loss. Solution of this equation gives space-time dependence of temperature. The finite 
element analysis starts with discretizing the film and the substrate into numbers of "8-node 
hexahedral" or "brick elements", the elements for the film sitting exactly over those for 
the substrate (see FIG Ej). Length and width of the film elements were same with those 
of the substrate element. Height of each substrate element was same with the thickness 
of the substrate while that of film varied depending on the position of the element due to 
roughness of the film surface and they were taken according to the topography of the pristine 



171 ] , the temperature distribution inside each 



film. Then following standard FEA procedure 
element was approximated as: 

T(x,y,z,t)=j2T?mf(x,y,z) (A.2) 

i=l 

where, Tf(t) are the temperatures at the nodes (nodal temperature) of the element at any 
simulation timestep, ipf(x,y,z) are some appropriate functions, known as shape function 
in FEA literature, with which temperature anywhere on and inside the element can be 
written. Using these approximation functions and following 'Galerkin weighted residual 
method' the equation IA.1I was re-written in a set of eight equations for each element 
where nodal temparatures were the unknown variables and the co-efficients were calculated 
using the shape functions and the thermo-physical parameters used in the equation lA.il The 
dimensions of the elements also enter into these co-efficients at this stage. Then, assembling 
all the elements from the film and the substrate, a set of equations was formed for all the 
nodes, solution of which gives the temperature at each node at any instant. Using the 
nodal temperature and the equation \A.2\ temperature anywhere inside the film and the 
substate could be calculated. Temperature at the next simulation time could be calculated 
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using forward time-step finite difference scheme. For this simulation we assumed that both 
the film and the substrate were thermally isotropic. We used p m = 1.049 x 10 5 kg.m~ 3 , 
s = 235.4 Jkg~ x K~ x and k = 430 Wm^R- 1 for silver film and p m = 2.330 x 10 4 kg.m' 3 , s = 
700 Jkg~ l K~ l and k = 150 Wm^ 1 K^ 1 for silicon substrate. At each simulation step, using 
the local temperature distribution obtained from FEA and calculating local stress differences, 
the atomic in-flux and out-flux for each film-element were calculated using the equation [31 
This gave the value of change in total number of atom inside the element. Multiplying this 
number by atomic volume, Q, change in volume of the element was calculated accordingly 
change in height also calculated (length and width of each element were kept constant 
throughout the simulation). 
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Figure Captions 

FIG [It (Color online) Schematic of experimantal set-up, pattern of the film is shown inset. 

FIG [2t Change in resistance (R) with time as the damage process progresses. The points 

marked by t Q , t\ etc are the stages for which SThM data are presented here. 

FIG 02 (Color online) (Row- wise) Topography, thermal map and typical temperature line 

profile (a) before current stressing, (b) after 146 hour and (c) after 163 hours for the part 

of film marked by the rectangle. All lengths are in fim and the temperature are in Kelvin 

scale. Dashed lines are to indicate the mean temperatures along the lines. 

FIG HI Evolution of (a) rms roughness (cr(h)) and (b) mean temperature (< T >) as a 

function of stressing time. Change in resistance also plotted to show the progress in the 

damage process. 

FIG [5t Change in (a) resistance (R), relative variances of temperature, (b) rms roughness 
of the film and local temperature gradient (VT) with time as the damage process progresses. 
FIG [6t (Color online) The film is modeled as a resistance network. For clean presentation 
only few resistance elements are shown. 

FIG E (Color online) Modeling of the film and substrate in numbers of finite elements. 
FIG [8} Evolution of (a) rms roughness (<j(h)) and (b) mean temperature (< T >) 
alongwith resistance as a function of stressing time with the progress of the damage process 
as obtained from the simulation. 

FIG [9t Change in (a) resistance (R), relative variances of temperature, (b) rms roughness 
of the film {cr(h)) and local temperature gradient (VT) with time as the damage process 
progresses as obtained from the simulation. 
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Pattern of the film 
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